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Abstract 

Measurement of Lepton-Flavor Violation (LFV) in the minimal SUSY Stan- 
dard Model (MSSM) at Large Hadron Collider (LHC) is studied based on a 
realistic simulation. We consider the LFV decay of the second-lightest neu- 



tralino, ~~ * ^' ~~ * H'Xii m the case where the flavor mixing exists in the 
right-handed sleptons. We scan the parameter space of the minimal super- 
gravity model (MSUGRA) and a more generic model in which we take the 
Higgsino mass \x as a free parameter. We find that the possibility of observing 
LFV at LHC is higher if \i is smaller than the MSUGRA prediction; the LFV 
search at LHC can cover the parameter range where the fi — > e^j decay can be 
suppressed by the cancellation among the diagrams for this case. 



1 Introduction 



The minimal supersymmetric (SUSY) standard model (MSSM) is one of 
the attractive extensions of the standard model (SM). Many experiments are 
searching for the possible evidence of the low-energy super symmetry. Among 
those, Lepton-Flavor Violation (LFV) processes may be considered as major 
discovery modes of supersymmetry; they do not exist in the SM or very small 
even if the small neutrino masses are introduced. 

In the MSSM the off-diagonal components of the slepton mass terms 
violate lepton-flavor conservation, and they are related to the origin of the 
SUSY breaking terms and interactions in physics beyond the MSSM. An 
approximated universality of the sfermion masses should be imposed at some 
energy scale so that the flavor-changing processes are suppressed below the 
experimental bounds. One of the candidates to realize the universality is 
the minimal supergravity (MSUGRA) model [Q]. In this model, the LFV 
slepton masses are induced by interactions above the GUT scale [f| or the 
right-handed neutrino scale ||. It is desirable to discover LFV in different 
processes so that we can reconstruct the off-diagonal slepton mass parameters 
which probe such interactions at the higher energy scale. 

Main constraints for the off-diagonal components of the slepton mass 
matrices come from the rare decay process searches at low energy, such as 
H — ► ej, /j,N — > eN, and r — ► y/y. The current experimental limits are 



Once the approximate universality of the scalar masses is imposed, the 
branching ratios are suppressed due to the GIM mechanism, and they can 
be smaller than these experimental bounds for reasonable SUSY parameter 
space. Therefore, those limits are not too serious at present. In future, some 
of the proposed experiments aim to reach to 



Br ((j, -> ery) < 1.2 x KT 11 [4], 
Br(fiN -»■ eN) < 6.1 x 1(T 13 §, 
Br{r -> p 7 ) < 1.0 x 1(T 6 @. 



(1) 



Br(n — > ej) 
Br(fxN -> eN) 




[7], 
1, 



2 



< io- 18 g, 



LFV may be searched for at future collider experiments through the pro- 
duction and decay of the slepton. While the rare lepton decay widths suffer 
the suppression of the order of (Am ; ~/m z ~) 2 , the production and decay pro- 
cesses of the slepton merely receive the suppression of the order of (Amj/Tj) 2 



12| . Here, Amj is the mass difference between the sleptons, and m z ~ and r ; ~ 
are the average mass and decay width of the sleptons, respectively. Thus, 
the future high-energy collider experiments could explore the region of the 
parameter space which may not be reached to by the rare decay searches. 
At future e + e~ linear collider experiments, the e production cross section 
could be very large if the bino-like neutralino mass Mi is relatively light, and 
then the e-\i{r) mixing may be discovered there ||12|| ||13| . Similarly, a muon 
collider has a potential to access to /x-r mixing |TJJ]. 

For the LFV search at Large Hadron Collider (LHC), the useful mode is 
the LFV decay of the second-lightest neutralino (x«D> X2 ~ * H ~ > Xi^'i where 
the slepton oscillation effect leads to the different flavors for the leptons 
in the final state (I and I'). This mode has an advantage to observe LFV 
compared with the direct Drell-Yang production of the sleptons, since x\ can 
be copiously produced through the cascade decay of the squarks and gluinos 
]15| . Typically 60% of the first- and second-generation left-handed squark 



decays into the wino-like neutralino and chargino, and in various models the 
right-handed slepton masses are predicted to be smaller than the second- 
lightest neutralino mass. By using the above mode, LFV in e-p, mixing has 



been investigated by Agashe and Graesser [15|, in which they chose the point 



5 of ATLAS TDR study pi in the MSUGRA model. The f-fi mixing has 



been recently studied by Hinchliffe and Paige fP7 |. 

In this paper, we estimate the reach of LFV at LHC based on a realistic 
simulation in the MSUGRA as well as a more generic model. We assume 
that LFV comes from mixing of the right-handed smuon and selectron. The 
signal of LFV \\ decay is the two opposite-sign leptons (e + /i~ or e~fi + ) in the 
final state. The distribution of the lepton-pair invariant mass (m efl ) in the 
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LFV final state has an edge whose value is known because it is the same as 
the edge of the lepton-flavor conserving opposite-sign modes (m ee and m M/1 ). 
Since the distribution of the background processes does not have an edge, we 
can obtain a sizable S/N ratio in the region near the edge of m eil . 

In the MSUGRA model, LHC can reach to M ~ 400 GeV where M is 
the common gaugino mass, when the mass difference and mixing angle of 
the right-handed selectron and smuon are Am =1.2 GeV and sin 29 = 0.5. 
Unfortunately, the broad parameter space in the MSUGRA, which can be 
probed by LHC, are already excluded by /i — > constraint. We point out 
that this is only for a case /i » M where /x is the Higgsino mass. The \x 
parameter may be close to the gaugino masses in the more generic model, 
and in the case the reduction of Br(fi — > ej) by cancellation among the 
diagrams is generic. We study a model of this kind (cMSSM) and we find 
that the experimental reach of LFV search at LHC extends up to M ~ 500 
GeV for the same slepton oscillation parameter. 

We organize this article as follows. In the next section we discuss prospect 
of the LFV searches at LHC in the MSUGRA and the cMSSM, and present 
sensitivity at LHC for generic parameters in Section 3. Section 4 is devoted 
to conclusions and discussion. In Appendix A we present formula we used in 
this paper, and in Appendix B the mass spectrum and the branching ratios 
for the SUSY particles are given for the sample parameter sets we adopt in 
this paper. 

2 Prospect of LFV searches at LHC in the 
MSUGRA and the more generic models 

In the MSUGRA model, the GUT scale scalar mass m and triliniear coupling 
A are universal among generations, and the gaugino masses in the MSSM are 
given by the SU(5) gaugino mass M. The universal scalar mass is predicted 
when the SUSY breaking sector couples to chiral multiplets equally. However, 
it is possible, and might be natural, to assume the different scalar masses 
between the matter fields and the Higgs fields in the GUT context. For 
example, in the most simple SO(10) model, we can set the following boundary 
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condition, 

m Q = m q R = m i L = m I R (= m ie) ; 

m Hl = m H2 (= m w ). (2) 

We call this choice as the cMSSM. The MSUGRA model predicts // > M. 
On the other hand, the cMSSM allows /i comparable to M. We will discuss 
the phenomenology of these models in the context of the LFV search at LHC. 



2.1 The MSUGRA model 

We summarize qualitative features of the MSUGRA predictions relevant to 
the LFV study at LHC. 

• fj, > M 

The Higgs masses at the GUT scale are common to the other scalar 
masses. The Higgsino mass parameter \i can be calculated as a func- 
tion of m, M, and tan/5 so as to reproduce the correct electroweak 
symmetry breaking. Especially for m t ~ 175 GeV, the dependence on 
m effectively disappears when m is of the order of M, which is some- 



times referred to the focus point [18] . In such fx/M is more or 



less fixed, and the LSP (xi)> the lighter chargino (xf), and the second- 
lightest neutralino are gaugino-like. The fi parameter is lighter 
than the gaugino masses only when m 3> M, while the decay of to 
I, which is relevant to our study, is closed in the case. Thus, one can 
safely assume x\ an d Xi are wino- or bino-like in our study. 

Light fx 

The lighter stau f\ is lighter than the other sleptons due to the left- 
right mixing proportional to tan /3 x fi. Note that the lower limit of the 
Higgs mass now strongly constrain small tan (3 cases |H| , and therefore 
the mixing is expected to be significant. The typical outcome is an 
increased branching ratio of x\ to f, and reduced branching fractions 
to the first- and the second-generation sleptons Brix® ~^ P)- 
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The fact that the MSUGRA model tends to reduce the branching ratio of 
X2 to the first- and second-generation sfermions directly limits potential to 
search for LFV at LHC This can be seen in Fig. [TJ, where we plot contours 
of some constant branching ratios of X2 f° r tan/? = 10 and fi = 1.5M 2 in 
the cMSSM, where M 2 is the gaugino mass of SU(2)x. The /i/M 2 ratio is 
roughly what the MSUGRA model predicts for moderate tan (3 and M ~ m. 
Due to the Higgs mass constraint, we do not plot the region where M < 300 
GeV. It is found from this figure that the branching ratios to cr or p,R are 
less than 6% for M > 420GeV. The reasons are the following. First, note 
that the decay into the Higgs boson opens for M > 325 GeV, and it quickly 
dominates over the other decay processes. Since x 2 * s wino-like, the decay 
X2 — * hXi 1S suppressed by only Nh oc 1//x, while x 2 ~^ * s suppressed by 
Nb oc l//i 2 . Here Nh and Nb is the Higgsino and the bino components of X21 
respectively. As \i is relatively high in the MSUGRA model, the branching 
ratio to h is significantly larger than that to Ir. Second, large /ztan/3 induces 
non-negligible f left-right mixing. Due to the left-right mixing, x\ decays into 
f\ through the tl component in fi, and the decay branching ratio dominates 
over the branching ratios to or jiR. Third, if m M, the decay into II 
is open, and this quickly dominates over the other decay modes. Thus, the 
region where the decay of x® to cr or p,R has a sizable branching ratio is 
limited. We will discuss the LFV signal in the next section. 

2.2 The cMSSM with fi ~ M 2 

We now turn into the case where universality of the sfermion masses does not 
hold for the Higgs sector as Eq. (g). We first argue that this model predicts 
different phenomenology for the LFV search. 

In Figs. and ^, we show the x\ branching ratios for [i = M2 and tan (3 = 
10, 20. Not only the decay to cr is kinematically open, the branching ratio 
is larger than that of the MSUGRA case. This is because the relatively large 
Higgsino-gaugino mixing makes the X2 mas s smaller than M 2 and the mass 
difference between x\ an d Xx becomes smaller. This makes the decay of X2 
to h close up to M ~ 420 GeV. At the same time, the decay into II is not 



6 



BR(X 2 -> X ) 1^=1-5 M 2 , tanp=10 




100 150 200 250 300 



m 16 [GeV] 

Figure 1: Contours of the constant branching ratios in the cMSSM model for tan (3 = 10 
and [i = 1.5M 2 . Solid lines are for Br(x2 -> hxx) = 0.3, 0.5, 0.7, 0.9 (from left to right), 
dotted lines for Br(yQ ~> ttj.) = 0.7, 0.5 0.3 (from bottom left), and thick solid lines for 
Br(x2 — * Gin) = 0.05, 0.03 (from bottom left). No solution above the thick line at the 
top left. 

open in the Figs. || and ||, because = m\ & + 0.8Mf — 0.2m| ^ M| > m| . 
The decay into T\ is also suppressed due to the reduced left-right mixing. 
Thus, in the cMSSM with fi ~ M, the prospect of finding LFV at LHC is 
considerably better compared with the MSUGRA model. 

Varying the /i parameter also non-trivially changes the fi — > cy branching 
ratio. In Fig. |], we show the contours of constant Br(fi — > ej) in m\§ and 
[ijM.2 plane for the GUT scale boundary condition given in Eq. (|2|). Here 
we take tan (3 = 10 and M = 300 GeV. We assume the only non-zero mixing 
mass term between cr and fiR exists at the GUT scale. We take the difference 
of scalar masses at the GUT scale and the mixing as Am = 1.2 GeV and 
sin 29 = 0.5. 

In the MSUGRA model, which corresponds to /VM2 ~ 1-5, Br(fi —>■ ej) is 
less than 10 -11 for itliq = m > 210 GeV and becomes minimum for m\§ ~ 300 
GeV. The suppression around mi 6 ~ 300 GeV is due to cancellation among 
diagrams that will be discussed soon. Although LFV in the sfermion decays 
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does not suffer from such a cancellation among diagrams, the decay of x° to 
I is closed for m\% = m ^ 210 GeV. This shows that it is difficult to observe 
LFV at LHC in the MSUGRA model. The constraint from Br([i — > ery) is 
strong to observe LFV at LHC in the MSUGRA model for generic slepton 
oscillation parameters as we see in the next section. 

The constraint weakens in the region where x® — > IrI is open, when \i 
is smaller than the MSUGRA prediction. The parameter dependence of the 
cancellation can be explained by the mass- insertion formula f2"0fl , which is 



expressed by the off-diagonal component ^i R f lR in the right-handed slepton 



mass matrix. It occurs among the 4 different amplitudes where the chirality 
flip occurs on either in the external or internal lines. Among them two dia- 
grams involving the left-right mixing or lepton-slepton-Higgsino vertex have a 
common overall factor m \ R p, R l m \ R x Mifi tan /3/m~ R with opposite-sign coef- 
ficients. When \i tan (3 is large and the absolute values of the two amplitudes 
are larger than the others, a nearly complete cancellation could occurs when 
fi/mj R ~ ///mi6 ~ 1.5, as in Fig. [|. When fi < 1.2M 2 , the cancellation 
occurs for m\e < 200GeV where the decay of \2 t° Ir is kinematically open. 

We have been discussing the case where LFV occurs due to the non-zero 
m^ fl . We will comment on the other cases in the last section. 

3 Potential of the LFV search at LHC in the 
cMSSM model 

In the limit of lepton-flavor conservation, the process X2 ~ * M ~~ > ^Xo a ^ LHC 



has been studied by many authors [16]. The distribution of the lepton-pair 



invariant mass (mu) in the final state is given by 



dTi-i+ j m n (0 < m u < m^ ax ) 



dmu 1 
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where the edge m™ ax is expressed by the slepton mass m~ { and the neutralino 



masses m^o as follows: 



ml rn 2 ~ 



(mr x ) 2 = m|(l-^)(l-^). (4) 

X2 m~ m~ 
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BR(% 2 l^=M 2 , tanp=10 




m.JGeV] 



Figure 2: Contours of constant branching ratios in the cMSSM model for tan/3 =10 
and [i — Mi. Solid lines are for Br(x2 ~ * ^Xi) = 0.5,0.7 (from left to right), dotted lines 
for Br(x2 — * tt®) = 0.6,0.2 (from bottom left), and thick solid lines for Br(x2 —> e &R) = 
0.2. The LSP is charged above the thick solid line at top left. 



BR(% 2 °^ X ) !^=M 2 , tanp=20 




100 150 200 250 300 



m 16 [GeV] 

Figure 3: Contours of constant branching ratios in the cMSSM model for tan/3 = 20 
and /j, — M 2 . Solid lines are for -Br(x° — > ^Xi) = 0.1,0.3,0.5 (from bottom to top), 
dotted lines for Br(x,2 ~~ * T ^i) = 0.9,0.8,0.6 (from bottom to top), a thick solid line for 
Br(\2 ~~ > e ^R) = 0.1. The LSP is charged above the thick solid line at left. 
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BR((X^ey) M=300GeV, tan|3=10 




100 150 200 250 300 
m.JGeV] 



Figure 4: Contours of constant Br(fi — ► ey) in mie and /V-M2 plane. A solid line 
corresponds to 10~ n , long-dashed lines 10~ 12 long and short dashed 10~ 13 , short-dashed 
10~ 14 . Here, M = 300GeV and tan/3 = 10. 

This decay process would be identified through the edge in Eq. (|3|). The main 
background comes from uncorrelated leptons from different squark or gluino 
decay chains. Fortunately, the background are estimated using the e ± /j^ 
distribution, and the distribution is smooth and decreases monotonically as 
m efJi increases. Thus, the background can be subtracted, and the subtracted 
distribution has a canonical structure as dT/dmu oc m« and terminates at 

iyy) mux 

m u ■ 

The signature of LFV on the process is the edge structure of e ± ^ dis- 
tribution on top of the accidental leptons from uncorrelated sources, since 
such an efficient subtraction method as above does not exist. The level of 
the signal and the background would be estimated if the production cross 
sections, acceptance of the signal and the background, and the background 
distribution could be estimated. This is of course possible by doing the MC 
simulation for each model parameters, however, we present semi-analytical 
approach in this section. In the following subsections, we will discuss the 
production cross section of the SUSY particles, level of the signal and the 
background, and the background distribution, and show the experimental 
reach of the LFV search at LHC at the end. 
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M(GeV) 


m 


m- g 


rrig L 


mi 


A) 


300 


100 


706 


633 


470 


B) 


350 


125 


818 


729 


549 


C) 


400 


150 


915 


824 


627 


D) 


500 


175 


1135 


1012 


781 



Table 1: Description of the set of parameters A)-D). We chose tan/3 = 10 
and A = in the MSUGRA model. The GUT scale gaugino mass M, the 
scalar mass m, the squark and gluino masses are given. 



(pb) 


er(total) 






o(qq) 


a{qq*) 


A) 


25.0 


3.03 


13.04 


4.70 


2.89 


B) 


11.2 


1.11 


5.67 


2.41 


1.41 


C) 


5.74 


0.45 


2.87 


1.44 


0.67 


D) 


1.60 


0.09 


0.76 


0.51 


0.19 



Table 2: Production cross sections for the parameters A)-D). 



3.1 Production cross section of the SUSY particles 

We start our discussion from an estimation of the squark and gluino produc- 
tion cross sections. The second-lightest neutralino x!> i s produced through 
the cascade decays of qi or g. The signal rate of the LFV decay x° ~^ ^Xi 
depends on the production cross section significantly, because it reduces very 
quickly with increase of the squark and gluino masses. 



We use the ISAJET version 7.51 |2T[ to estimate the production cross 
sections.[] We are interested in the region of parameter space where m < M so 
that the x° decay to / is open. In this range, we choose the four parameter sets 
A)-D) in Table 1 as samples and derive a fitting functions for the production 
cross sections in Table 2.0 The fitting functions will be used later to estimate 
the number of the signal and background e /i T events. 

Since the squark and gluino masses are quite close to each other for those 
points, we fit the production cross sections by the following simple functions 

lr The parton distribution is given by CTEQ3L. 

2 We omit the cross sections of the stop- and sbottom-pair production since they give 
negligible contributions to the jie background. 
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(j(total) 


<r{$9) 


<r(9Q)+<r(9Q*) 


a(qq) 


a(qq*) 


ax 


3.34 


0.23 


1.74 


0.92 


0.39 


a 2 


5.78 


7.407 


5.780 


4.677 


5.733 



Table 3: Results of the fit of the production cross sections to the functions 
given in Eq. (0). 

of nig, 

a(qoig) = a 1 (m~g/TeV)~ a2 (pb). (5) 

The production cross sections are fitted very nicely by the parameters a\ and 
a,2, listed in the Table 3. When the scalar masses are changed in a relevant 
region of the parameter space, the cross sections are changed within only 
20%, which is within QCD uncertainties. 

Note that the scaling parameter a 2 depends on the production modes 
substantially. Processes involving initial state gluon(s) such as cr(gq), <y(q*g) 
reduces quickly compared with a(qq) production when gluino and squark 
masses become heavier. This comes from a fact that the parton distributions 
of the gluon and sea quarks are generally softer than valence quarks. 

The u/d and qi/qn ratios affect the level of fie background. The u^d^, 
ulu* l , and d^d\ production are likely to contribute to the signal compared 
with the other squark productions, because the up- type squark ul would 
decay into xf producing l + while d L producing l~ , and q^ decay mostly 
into Xi when Xi, Xij xX are gaugino-like. Also, note that the production is 
dominated by valence u and d components of parton distribution function 
as in Table 2. Thus, we adopt a simple approximation that u : d = 2 : 1, 
u* : d* — 1 : 1, and '■ qn = 1 : 1 in any squark and/or anti-squark 
production cross sections.^ 



In the MC simulation, the ratio of production cross sections, a(uu) : a(ud) : <r{dd), 
is 4.4:3.9:1 for point A) and and 7:5.3:1 for point D). For qq* production cross sections, 
the ratio a(uu*) : a(ud*) : <j(du*) : a(dd*) is 5.7:4.2:1:3.9 for point A). And, a(qLqL) : 
a{q L q R ) : a(q R q R )=l:lA:1.2 and cr(q L q* L ) : a(q L q R ) ■ cr(q* L q L ) : cr(qRqR)= 1:0.93:0.90:1.15 
for point A). 
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3.2 Level of the signal and the background and the 
acceptance 

Next, we will estimate the level of the signal and the background, and eval- 
uate the acceptance. For this purpose, we use the MC data for the following 
parameters^; 

point I) (1 = 497.87GeV, tan/3 = 2.1, M = 300 GeV, and m = 100 GeV in 



the MSUGRA \22\ 



point II) n = 199.85GeV tan/? = 10, M = 250 GeV, and m w = 90 GeV in 



thecMSSM 23 



To estimate level of the background, we have to calculate branching ratios 
Br(g -> l±X), Br(g -> l + l'~X), Br(q -> l±X), and Br(q -> l+l'-X). We 
find that g has a large branching ratio into multiple leptons. For example, 
g tt and g — > &i& followed by bi — > t\ + or b\ — > tW 7- have fractions 15%, 
6%, and 4%, respectively, for point I). This results in a high probability to 
get multiple leptons in the final states; Br(g — > (t,b 1 ,b 2 ) — > l) — 4% and 
Br(g — >■ (t, & 1; 6 2 ) — >■ i'/ or lr(— > /')) = 2%. Thus, the ^ production may be a 
significant source for the background for the LFV search compared with the 
squark-pair production processes, since the squark-pair production processes 
are required to have two cascade decays involving a lepton each so that they 
contribute to the background. 

The chargino production is also a significant source of the background. A 
chargino may decay into W, u, and f 1; followed by decay into leptons in the 
final state.Q A produced tau lepton may further decays into e or fi, whose 
branching fraction is about 35%. The second-lightest neutralino %2 may also 
decay into rf\. While the momentum of e or fi from the tau-lepton decay 
tends to be low, the acceptance of such leptons would not be negligible for 
the adopted lepton px cut as low as 10 GeV. For the case where \i ~ M 2 , one 
should also take care of the decay of q into xlu)! xf- They are also calculated 
and added in our background estimation. 

4 The mass parameters and relevant branching ratios are given in the Appendix B, 
which are calculated by the IS A JET. 

5 For example, 89% of chargino decays into fi for point I). 
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The total e ± fi T X and l + l~X events before cuts are shown in Table 4. The 
numbers will be compared with the MC results in order to derive the accep- 
tance. They are obtained by calculating the gluino and squark branching 
ratios into the lepton(s) using the major branching modes mentioned above, 
and multiplying our fitted production cross sections described previously to 
the branching ratios. This semi-analytical calculation is checked by inde- 
pendent toy MC simulations which includes all decay steps. We include the 
contribution from qq, qq* and qg production. Nw is the number of the e ± /i =F 
events coming from primary leptons from W or I decay. Ni T is the number 
of the events with one primary leptons and one lepton from r — > e or //. 
N TT is the number of e ± / u =F events from the leptonic decay of two r's; here 
we omit the contribution from x°i r+T ~Xii because the mm distribution 
are substantially softer than expected LFV signal x° — > l' 'I ~^ ^Xi- Note 
that our calculation does not include the probability that 3/, or 4/ events are 
accepted as two lepton events etc. 

Point II) is substantially lepton rich compared with point I) due to the 
enhanced branching ratio into sleptons. This is typical for parameters with 
H ~ M 2 with light sleptons as shown in Fig. |2|. 

Now we can estimate the acceptance of events involving , e + /i~, or 
e~/i + . Here we use MC data produced and simulated by the IS A JET and 
the ATLFAST. 2 x 10 6 (1 x 10 7 ) events are generated for point I) (point II)) 
corresponding to 95(196)fb _1 of the integrated luminosity. In the simulation, 



we adopted the cuts given in |24| 



• E^ iss > max (lOOGeV, 0.2M eff ) , 

• Pj x > 100 GeV, Pf 2thth > 50GeV, 

• M eS > 400 GeV, 

• a pair of isolated opposite sign leptons with Pf > 10 GeV. 

Those cuts are chosen to reduce backgrounds from QCD processes to a neg- 
ligible level. After the cuts, the e ± if iL events are reduced to 3609 events for 
point I), and 42721 events for point II). (See A r (e//)(MC) in Table 4.) By 
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point 


7VM (MC) 


N(l'l) 


N(lr) 


N(rr) 


point I) 


3609 


1.49 x 10 4 


3.61 x 10 3 


3.2 x 10 2 


point II) 


42721 


3.1 x 10 4 


10.8 x 10 4 


5.2 x 10 4 


point I) b veto 




1.58 x 10 3 


5.4 x 10 2 


0.67 x 10 2 


point II) b veto 




3.78 x 10 3 


3.85 x 10 4 


2.20 x 10 4 



Table 4: Number of events with e ± /x =F before cuts from different decay processes N(l'l), 
N(It) and N(tt) when lepton flavor is conserved. (See text.) The numbers are for lepton- 
flavor conservation. JV(efi)(MC) is the numbers of accepted events with e ± /i T satisfying 
0GeV< m eM < 200GeV, obtained by the 1SAJET+ATLFAST simulation. 

comparing N(l'l), N(It), N(tt) with JV(e/i)(MC), we estimate the accep- 
tance of uncorrelated e/j, events at the ATLAS detector; 19% for point I), 
and 22% for point II). 

The acceptance of LFV signal may be also estimated by looking at the 
the acceptance of opposite-sign same-flavor leptons in the same simulation. 
Note that the LFV signal from the x® decay has exactly the same kinematics 
to the lepton-flavor-conserving decay. We list in Table 5 the result of the MC 
simulation of the same samples as that of Table 4 and the estimated 
production from the x\ decay before the cut. The acceptance of the process 
is 28%. This is higher than that of e ± fi T events under the same cut. To be 
conservative, we adopt constant 25% acceptance for both the signal and the 
background. 

To suppress the background furthermore, the 6-jet veto may be efficient 
while we will not use it in this article. In the Tables 4 and 5, we show numbers 
of events which do not involve the third-generation squarks in the cascade 
decays. The events involving t or b could be removed by the 6-jet veto. If 
the efficiency of the veto is ideal, N(l'l) is reduced by an order as in Table 4, 
while the number of the correlated leptons from x\ I decays does not 
change in Table 5 by the 6-jet veto. However, the rejection factor depends on 
the decay patterns. For example, xt ~^ t dominates the chargino decay for 
point II), and the decay — > xt becomes the efficient source of r lepton. 
This means that N(lr) is not reduced so much by the 6-jet vet in point II). 
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point 


/ + / events before cut 


accepted l + l by MC 


point I) (951b- 1 ) 
point II) (196 fb- 1 ) 


8.72 x 10 4 
8.97 x 10 5 


2.55 x 10 4 
2.55 x 10 5 


point I) b veto 
point II) b veto 


6.51 x 10 4 
6.96 x 10 5 





Table 5: The number of \2 ^ l + l Xi when LFV is absent. The accepted l + l events 
in MC simulations are also given in the table. 

3.3 Background distribution 

The signal distribution from decay increases with mm as in Eq. ([3]) and 
the distribution has an edge determined by the neutralino and the sfermion 
masses. On the other hand, backgrounds come from t, W, and xf decays and 
do not have the edge structure. They reduce rather quickly as increases. 
Therefore it is better to use the data near m™ ax so that S/N ratio maximizes. 
We should note that the position of the mm edge is known precisely from the 
same-flavor opposite-sign mu distribution. 

To estimate the number of the background events near the edge, we again 
use the MC data and fit it to the following fitting function, 

= k{mg) exp ( — — x mu) . (6) 
dm u V m g ) 

The background distributions and the fitting curves are shown in Figs. || and 

|. The data between 40 to 200 GeV and 100 to 200 GeV are used for the 

fit. The best fit is obtained c = 10.4 for point I) and c = 13.7 for point II) 

for 30 GeV < m eil < 200 GeV. The fitting function reproduces the large mm 

region reasonably well, but it fails significantly in the small mm region. It is 

natural that the distribution has a certain peak, which must be proportional 

to a typical momentum of leptons of the uncorrelated production process, 

such as the half of W boson mass. The distribution beyond this peak must 

be more sensitive to typical momentums of W, t, or xti which may depend 

on the gluino and squark masses. This is the reason we choose Eq. as 

a fitting function. The average value for c of those two points, c = 12.1, is 

used for our background estimation. For the plot, k(M) may be fixed so that 

overall normalization agrees for the region used for the fit. 
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Figure 5: The m efi distribution for point I). The data corresponds to integrated lu- 
minosity of 95 fb _1 and standard cuts are applied (see text). The histogram shows the 
distribution without LFV, while bars are number of events and the error with fl-e mixing. 
In the plot, 1/30 of x° ~~ > l"U ~~ * Xi'' decay chain is assumed to go to the e/Lt channel. 
Two curves are fits to the background distribution in the region ma = 40-200 GeV(solid) 
and mu = 100-200(dashed then solid). We use c = 12.1. 

mixed case 




50 100 150 200 

mJGeV] 

Figure 6: Same as Fig. |5], but for point II). The integrated luminosity is 196 
fb -1 , and c = 13.7. 
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3.4 Significance of LFV at LHC 

Having gone through all estimation needed, we now calculate significance 
of the LFV signal at LHC. Numbers of the signal (N sl9 ) and background 
(N b9 ) are estimated by (the fitted cross sections) x (the branching ratios) 
x (the overall acceptance 25%) x (the integrated luminosity). We assume 
the signal distribution in Eq. (|3|) and background distribution in Eq. (|J). 
We determine the overall normalization factor of Eq. (|6|) so that number of 
background above mm > 20 GeV agrees with the estimation. 

We define the A% 2 using the estimated signal and background events 
between 

max ^30GeV, < m u < m^ ax , (7) 

and calculates 



AX 2 = E 



iVf g + ;vt g 



for the bin size 2n GeV, where integer n is determined so that max(Af- lg ) > 
10. Eq. (|) expresses the statistical significance of the signal after the sub- 
traction of expected background.^ In the experimental situation, one may 
determine the background distribution from the real data, when the events 
above the m™ ax may be used to make a simple extrapolation as suggested in 
Eq. (D.Q 

In Fig. [7], we show the contours of constant A% 2 in mi6 and M plane. 
Ax 2 = 25 ~ 5cr contours correspond to 70 signal e ± n T events in the signal 
region in Eq. (0). Here, the integrated luminosity is 100fb _1 . The SUSY 
background is roughly of the same order as the signal. In Ref. PS|| , 120 total 
SM background events are expected for 30 fb _1 under the cuts 1) E^ iss > 300 

6 We assume no error for the background shape. 

7 Alternatively, one can constrain the MSSM parameters as model-independent as 
possible, so that the background distribution can be determined model-independently. 
Note that the nature of the third-generation sfermions is important since the substantial 
fraction of the e ± ^t T events might come from g — > tit or bb followed by their cascade 
decays to W etc. The sbottom mass may be reconstructed from the distribution of the 
events with b jet M&. Attempts to reconstruct stop decays may be found in flH, M 
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GeV, 2) pf > 10 GeV, and 3) two jets with pj et > 150 GeV. The level of 
the background in the signal region is of the same order as that of the SUSY 
background. The significance beyond this 5a contour is therefore the subject 
of more careful MC simulations both for SUSY and SM backgrounds. 

When [i ~ Mi, the parameter space covered by LHC extends, compared 
with the MSUGRA model (/i ~ 1.5M 2 ), due to the enhanced branching rates 
of X2 to Ir as we see in the previous section. We can see another qualitative 
difference in mis -C M region. For the MSUGRA case, the small mis region 
cannot be reached because x° ~~ > II dominates. The search region is extended 
to this region for the cMSSM with fi ~ M because m^o < M 2 and X2 could 
not decay into II for \i ~ M. 

We also estimate the LHC reach for generic oscillation parameters. In 
Fig. |8|, we plot the 5a contour and the line of Br(fi — > ej)—1.2 x 10~ n , 
1.0 x 10~ 12 , and 1.0 x 10 -14 in the parameter space of sin 28 and Am at 
the GUT scale. We use the the MSUGRA model with tan/5 = 10, A = 0, 
m = 100 GeV, and M = 300 GeV. We can see that the most part of the 
parameter range where LFV can be observed at LHC is already excluded 
by the Br(fi — ► e^) constraint, and the remained region will be covered by 
next-generation experiments. The situation changes for the cMSSM case with 
fi = Mi- The corresponding figure is shown in Fig. || Because of the change 
of the decay kinematics, the wider range of the parameter space is covered by 
the LHC than that in the MSUGRA case. Br(/j, — > ej), in contrast, becomes 
small due to the cancellation between the diagrams. It follows that LHC 
might be more advantageous to observe LFV than the \i — > decay search, 
especially in the cMSSM. 

4 Conclusions and Discussion 

In this paper we investigate the potential of LHC to find LFV in \\ ~^ 
I" I'Xv Here we studied it in a general model where the \x parameter is 
independent of the gaugino mass M by allowing the non-universal GUT scale 
Higgs masses. An approximated universality of squark and slepton masses is 
imposed. We find LHC would be able to find the LFV mixing between the 
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1 — r 




m.JGeV] 



Figure 7: y A% 2 = 5 contours for the LFV discovery. The thick solid line is for 
/i = 1.5M 2 and tan (3 = 10 in the cMSSM, the thick dashed line for fj, = M 2 and tan/3 = 20, 
and the solid line for fi = M 2 and tan/3 = 10. We fix the en-fin mixing angle 8 as 
sin 28 = 0.5 and the slepton mass difference Am = 1.2 GcV at the GUT scale. 



first and second generation in the right-handed slepton masses, mg R ji R while 
the £>r(/i — > ej) is undetectably small. 

LFV in left-handed slepton mass matrix as Ta\ L r L might be more moti- 
vated when the data from atmospheric neutrino is considered [14] [27]. We 
note that the cancellation among the LFV diagrams is unlikely when only 
the left-handed slepton mass is the source of LFV. We show the relation 
between Br(r — ► yrf) (the dashed line) and Br(x2 —> Xi r A i ) normalized by 
Br(x2 ~^ Xi rr ) (the solid line) in Fig. [T^ when TT& ^ is the unique source 
of LFV. Hinchliffe and Paige stated that it is possible to observe LFV for 
Brfy® ~^ Xi T f J ') ~ 0-01 i n this figure. We can see in Fig. 10 that the reach of 
LHC corresponds to Br{r — > ^7) ~ 10 -6 , which is also within the range of 
the r — > /X7 search at the KEKB experiment |TT . 



Finally when several off-diagonal scalar masses m~ ;; ~ are non-zero, Br (I — > 
l >r y) could show the very complicated structure, therefore negative results in 
the rare decay search at low energy do not necessary constrain the processes 
involving I decays at the high energy colliders. Especially, the LFV processes, 
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sin26 



Figure 8: The LHC reach and the line of the constant Br{pb — > ej) in the MSUGRA 
model are shown. Here, tan 0=W,A = O,m= 100 GeV, and M = 300 GeV. 




sin26 



Figure 9: The LHC reach and the line of the constant Br(p, — > ej) in the cMSSM are 
shown. Here, pt = M 2 , tan/3 = 10, A = 0, m 16 = 100 GeV, and M = 300 GeV. 
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10 




Figure 10: -Br(x° — ► Xi T M) normalized by Br(x° ~^ Xi TT ) (the solid line) and Br(r — > 
/ry) (the dashed line) as functions of m~ ifi /m 2 . Here, we take m = 130GcV, M = 
250GeV, A a = 300GeV and tan/3 = 10. 



/i — > cy and /xiV — ► eA^, may be more sensitive to all LFV slepton masses, 
compared with those involving the third- generation sleptons, such as r — > /ry. 
The size of m? f and m?- is less constrained compared with m?-, and then the 
LFV through f could overcome direct jx-e mixing. For example, even when 
only the left-handed sleptons have the LFV masses, there is a cancellation 
in fj, — > e7 among the diagrams in some specific parameter space 
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Appendix A:Formula for LFV processes 



vv Jq 2 >0 



First, we present our formula for the LFV decay of the second-lightest neu- 
tralino, ~^ —*)l + X°i an d -^)l'~Xi- The invariant mass (m v i) 

distribution of the leptons is given as 

dT 

(r&tf)Br$$tf) + r^tf)B4% 2 )) ■ (A.i) 

Here T x ^ Y (q 2 ), Br x ^ Y (q 2 ), and p(q 2 ) are 

, (A.2) 

m hJ 

2 2 / 772 - \ ^ 

Br ^ (9 } " lev \ R ^ R *y + L ^ L ^ / pf]77 I 1^ J (A } 



rW — ^2 f p(n) pW* i 7-H r( n )*\™ | 

i XYW J - \ U 12X U 12Y +%%| m ^ I 



(m{rV)^ (° ^ m ^ ^ (™r(? 2 )) 2 ) 



pw j I o > ( m r(? 2 )) 2 ) • 1 • j 

where [mr] xy = {m lx T lx + m- w Y- w ) /2, and 

KT% 2 )) 2 = m|o(l--^)(l --*-). (A.5) 
The interaction Lagrangian of fermion-sfermion-neutralino is written as 

Ant = -92¥a{RT2xPr + L^xPlWx + h-c, (A.6) 
and the coefficients are 

L iAX = —E;{[-{0 N )a2 ~ (0 N ) A1 t W ]U X ,i + — ^—^{0 N ) A3 U X ,i+3}, 

y Z ifly/ COS p 

Rflx = 4={ mi J O N )AsU x ,i + 2(0 N ) A1 t w U X:i+3 }. (A.7) 

V 2 771 w COS P 

The function of slepton momentum A XY (q 2 ) is 

i 5(q 2 — m~ ) + 5(o 2 — m~ ) 

A ^ = ^ " (A ' 8) 
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with xxy = {tti~ — m f y )/2[^r]xy- 

We can simplify above formula assuming two-flavor mixing of the right- 
handed slepton and the almost degenerate masses, 

dT 2T Br 



x l'l 



dm ll 



(m 



sm z 26 



(A.9) 



for < mfq < (m™ ax (m~)) 2 . Here, 



9 2 y 



16tt Xa 



m^o 1 — 



xi. 



Br 



9v 



2 

I 



2 \ 2 



^7r [mf] 



(A.10) 
(A.11) 



and m~ and sin 28 are the average mass and the mixing angle for the sleptons. 

Next, we present formula for the LFV lepton decays fi — > ey or r — > /ry. 
Those rates are also written in similar Lagrangian though, this time, contri- 
butions from chargino loops are also important. The Lagrangian involving 
chargino-slepton-lepton is given as 

C = -92l{L^ x P R + Rfl x P L )xA^x (A.12) 
where x^(A = 1, 2) is a chargino mass eigenstate. The coefficients are 



L 
R 



(c) 
iAX 

(c) 
iAX 



(0 R )AlUx,i, 

\f2m-w cos (5 



{Ol)a2U x , 



Then 



where 



nh - hi) 



— mf (\A L \ 2 + \A 

167T ZjVI ' ' 



R\2\ 



A R = A^ R + A^ R A L = A^ L + A^ L . 



(A.13) 



(A.14) 



(A.15) 



The coefficients in above equations are given as |29 
1 ] 



A {n)R 



32vr 2 m? 



it; AY-tt, 



' lA ' lx 6{l-x AX ) 4 
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A (c)l 



x(l - 6x AX + 3x 2 AX + 2x AX - 6x 2 AX InxAx] 



„(n) r(n)* m xj 



m z . (1 - re ax 



AX 



2x A x In xax) 



327r 2 m 



2 



p(c) p(c)* 1 

^AX^AX 6(1 _ Xax)4 



x (2 + 3a; ax - 5x AX + x AX + 6xax In a; ax) 

-(-3 + 4xax - ^ax -21iixax) 



, , , v m-+ 

i d( c ) r(c)* Xa 

^^iAX^jAX 



mij (1 - XAXy 



L<->.R- 



(A.16) 



B. Sample points 

In this paper, we used MC simulation data for two sample points to estimate 
the event distribution and the acceptance. We summarize the masses of 
SUSY particles and decay branching ratios here because they depends on 
choice of gauge couplings and so on. 

First, we list mass parameters and relevant SUSY-particle masses in GeV 
for the point I), which is studied in this paper. The IS A JET is used to 
generate this spectrum. 



m 


100.0 


tan/3 


2.1 










M 


300 




497.87 


Mi 


126.23 


M 2 


252.36 




157.21 




238.78 


rriu L 


654.11 


m d, 


657.18 


VTLy 


230.19 




121.52 




630.95 




628.38 


A2 


233.02 




499.18 




232.03 


m xt 


520.11 




523.23 




230.14 




459.55 


m h 


670.68 




156.81 


rrif 2 


238.92 




600.30 


m h 


628.84 


m h 


94.33 


mp 


606.79 




611.67 


m 9 


732.94 



For the parameter gluino could decay into squarks, and especially, it has 
enhanced branching ratios to the third-generation SUSY particles. The first- 
generation SUSY particles would be also generated and decay into chargino 
or neutralinos to produce leptons. The major branching ratios in (%) are 
following; 
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o — > 61 b 


15 


— > 696 


9.7 


Q — > tit 


15 


Mr — » v2m 

i-v A, 2 


33 


V, T — > Y~!~d 
U 'L Al 


65 


u.^ A2" 


31 


d r — > 11 

U L Al u 


64 


^1 Al 1, 


23 


h -> $t 


12 


*! -> *+& 


64 


61 - 


28 




41 


61 -> t x w- 


27 




67 


b 2 -> Xf * 


9.6 




15 


X° -»■ e « e 


9.2 


X2 -»■ 


9.2 


x° -»■ f i r 


12.1 




89.7 




1.4 




0.5 


Xi~ -»• Tl^T 


6.6 







For point II), the input mass parameters and the SUSY particle masses 
are following; 



77li6 


90.0 


tan/5 


10 










M 


250 


n 


199.85 


Mi 


103.9 


M 2 


208.75 


m e R 


139.3 


m- eL 


206.09 




556.07 


m dr. 


561.67 


m 9 


190.28 




93.18 




534.36 


m d H 


533.21 


777^0 


155.13 




208.74 


m x\ 


148.44 


m xt 


272.52 


m x° 4 


273.8 




188.67 




374.43 




563.81 


m fl 


132.56 




206.06 


™>b, 


498.27 




531.2 


m h 


112.59 


ITLp 


436.98 


m H 


437.63 


rrig 


624.36 



Relevant branching ratios in % for the sample parameter are 
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q — > bib 


17.4 


5 — ► 626 


10.1 


g — > t x t 


13.2 


— > x 2 M 


20.7 




0.4 


-> X4 M 


12.2 


-> X2 rf 


21.4 


<4 -> xid 


14.5 


- * J&d 


0.7 


dh — >■ xV 


15.4 


C? r — > V^M 

i/ Az 


36.4 


*i -»■ *?* 

1 A, 1 


9.5 




10.1 


-L /VI 


74.3 


A2 


6.1 


L _> y+£ 
u l Al ^ 


29.9 




37.8 




9.8 




10.0 


&2 -> Xl"* 


14.1 


^2 -> X2 * 


35.5 


&2 -> Xo & 


4.9 


6 2 -> tiW" 


10.5 


X° -»■ g Le 


5.5 


X4 g /?e 


1.1 


X 2 -»■ e fi e 


23.6 


X2 -> i>Le 


9.7 


Xi~ -> 


89.1 


-»■ xt e 


40.1 


e L -> 


38.2 


e L -> x?e 


19.8 


e/? -> X?e 


100 
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